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ABSTRACT 

Using the Submillimeter Array we report the discovery of a compact low mass 
bipolar molecular outflow from L1014-IRS and confirm its association with the L1014 
dense core at 200 pc. Consequently, L1014-IRS is the lowest luminosity (L ~ O.O9L ) 
and perhaps the lowest mass source known to be driving a bipolar molecular outflow, 
which is one of the smallest known in size (~500 AU), mass (< 10" 4 M ), and ener- 
getics (e.g., force < 10" 7 M Q kms" 1 yr" 1 ). 

Subject headings: ISM: individual (L1014, L1014-IRS) - ISM: jets and outflows - 
stars: formation - stars: low-mass, brown-dwarfs - techniques: interferometric 



1. Introduction 

Spitzer Space Telescope observations of the Bok globule L1014, previously classified as star- 
less (Parker 1988), discovered an infrared point source with protostellar colors, L1014-IRS, located 
near to the globule's dust emission peak (Young et al. 2004). The infrared colors of L1014-IRS are 
consistent with (i) a protostar or proto-brown-dwarf in L1014 at ~200 pc (Leung, Kutner & Mead 
1982), as assumed by Young et al. (2004), or (ii) an intermediate mass T Tauri star (L ~ 16 L ) 
associated with a background cloud in the Perseus arm at 2.6 kpc. 

Assuming that the distance to L1014-IRS is 200 pc, its star plus disk luminosity, assuming 
isotropic emission, is ~0.09 L Q , and the best fit to its SED is with a protostar luminosity of 0.025 
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L Q . If L1014-IRS has the age of a typical Class I protostar of ~ 10 5 yr then it has a substellar 
mass of only 20-45 M Jup (Huard et al. 2005). However, no spectroscopic observations exist to 
determine its mass and accretion rate, and so it is unclear whether L1014-IRS is a young protostar 
still acquiring a significant fraction of its final mass, or whether it is destined to remain substellar. 

More recent observations strongly support the conclusion that L1014-IRS is embedded within 
L1014 (Crapsi et al. 2005; Huard et al. 2005), however the evidence is still circumstantial. It has 
not been shown that L1014-IRS has the same local velocity as L1014, nor has an outflow been 
detected, despite recent searches using CO with single dish telescopes (Visser et al. 2001; Craspi 
et al. 2005). In this paper we present Submillimeter Array (SMA) observations of CO 7 = 2^1 
toward L1014-IRS, to search for a compact outflow and to resolve the distance ambiguity. 



2. Observations 

Observations of L1014 near 230 GHz were obtained with the SMA 1 (Ho et al. 2004) on 2004 
August 15 and September 5. Zenith opacities at 225 GHz were typically in the range 0.1-0.2. 
The observations utilized both 2 GHz wide receiver sidebands, separated by 10 GHz. The SMA 
correlator was configured with high spectral resolution bands of 512 channels over 104 MHz for 
the 12 CO, 13 CO and C 18 7 = 2^1 lines, providing a channel spacing of 0.26 km s" 1 , with a lower 
resolution of 3.25 MHz/channel over the remainder of each sideband. Observations of L1014 were 
interleaved with the quasars BL Lac and J20 13+370 for gain calibration. The data were edited and 
calibrated using the MIR software package adapted for the SMA. Saturn and Uranus were used 
for passband and flux calibration, respectively. The flux of BL Lac measured on the two days 
agrees to within 20%, which we take to be the uncertainty in the absolute flux scale. Mapping 
was performed with the MIRIAD package, resulting in an angular resolution of l'.'2 x l."0 (natural 
weighting). The rms sensitivity was ~2 mJy for the continuum, using both sidebands (avoiding 
the band containing the CO line), and ~0. 1 Jy beam" 1 per channel for the line data. The primary 
FWHM beam of the SMA is ~55" at these frequencies. 



3. L1014-IRS: A Protostar at 2.6 kpc? 

The systemic velocity of L1014 is 4.2 kms" 1 (Crapsi et al. 2005). Behind L1014 lies an- 
other cloud with a velocity near -40 km s" 1 associated with the Perseus arm at 2.6 kpc (Young 
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et al. 2004; Crapsi et al. 2005). Emission from the 7 = 1 -> lines of C 18 and N 2 H + , and the 
7 = 2^1 line of CS, tracing gas density 10 3 " 5 cm" 3 , was not detected at these velocities in a ~50" 
beam by Crapsi et al. (2005). No CO emission is detected in the SMA data over the velocity 
range [-35,-45] km s" 1 , implying that no compact CO emission nor outflow is associated with this 
velocity component to the level of our sensitivity. 

A tentative detection of continuum emission of 7 ± 3 mJy was obtained at the position of 
L1014-IRS (offset by 4"8, 0."1 from the phase center) with observations made in the compact 
configuration, but was not confirmed by the extended array observations. By combining the two 
data sets and only selecting baseline lengths < 50 kA a similar tentative detection (5 ± 2 mJy) 
is obtained. Young et al. (2004) predict a flux of ~l-2 mJy at 230 GHz for the central object 
(star+disk) of L1014, which is consistent with our results. Recent BIMA imaging at ~95 GHz with 
a 15-20" synthesised beam shows weak slightly extended emission (Lai et al. in preparation), but 
no compact component. Young et al. (2004) predict that the emission from a background protostar 
(L ~ 16L ) would be ~ 0.6 Jy, which is clearly ruled out by the SMA observations, and together 
with the lack of line emission at -40 kms" 1 (§3), rules out the possibility that L1014-IRS is a 
distant young star. 



4. CO Outflow Associated with L1014-IRS 

CO 7 = 2^1 (hereafter referred to as CO) emission was detected at velocities blue-shifted 
and red-shifted with respect to the systemic velocity of L1014 of 4.2 kms" 1 . No emission was 
detected at the systemic velocity. Figure 1 presents an overview of the CO emission integrated 
over the blue- and red-shifted velocities, and Fig. 2 presents a position- velocity (P-V) diagram 
along the main axis of the compact outflow (indicated on Fig. 1(b)). The CO emission is clearly 
offset from the position of L1014-IRS, in a pattern typical of bipolar molecular outflows, with 
the blue-shifted emission aligned with the near-infrared scattered light nebula (Huard et al. 2005). 
Two representations of the data are shown. In Fig. 1(a) the data have been tapered with a 3" 
Gaussian in order to improve the sensitivity to the fainter emission seen ~ 20" north and south 
of the main outflow. The red-shifted emission to the south is clearly visible even at the highest 
angular resolution. These faint extended emission features are approximately equidistant from 
the position of L1014-IRS and may be associated with the outflow, perhaps indicating episodic 
emission. Fig. 1(b) and Fig. 2 show that the inner bipolar emission is very compact: the peaks 
have projected separation of ~200 AU centered on L1014-IRS, the lowest contours are separated 
by at most ~700 AU, and the size of each lobe is only ~540 AU (0.0026 pc) in length. 

No emission is detected at the cloud systemic velocity, and the peak emission in both the blue 
and red lobes are offset by similar amounts from the systemic velocity (Fig. 2). The red-shifted 
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emission is clearly brighter than the blue-shifted emission by about a factor of two. The brightness 
asymmetry might arise because L1014-IRS is offset from the density peak, with the red-shifted 
outflow propagating into denser material resulting in a larger swept up mass, compared to the 
blue-shifted emission. The core's column density peak lies ~ 10" to the south (red-shifted side) of 
L1014-IRS (Huard et al. 2005), supporting this conjecture. Alternatively, the difference could be 
due to radiative transfer effects implying a moderate CO opacity (Bally & Lada 1983). 

4.1. Molecular Outflow Mass, Kinematics and Dynamics 

In order to derive the outflow properties we only consider the compact outflow shown in 
Fig. 1(b). The nature of the "extended outflow" component is unclear, and requires confirmation, 
perhaps by combining the SMA data with single-dish data to recover extended missing flux. Ob- 
servations with the IRAM 30-m telescope of the L1014 cloud in CO 7 = 2^1 have recently been 
made by us. A full analysis of these data are beyond the scope of this paper. However, we can 
investigate qualitatively the missing flux question. Figure 3 presents a comparison of the SMA and 
30-m data. The SMA data have been convolved with a 4" beam (i.e., the approximate area of the 
outflow), and scaled by the ratio of the beam areas, to approximate the outflow emission detectable 
by the 30-m in equivalent brightness temperature units. Figures 3(a) and (b) show that no emission 
and so no missing flux is present at velocities greater than that detected by the SMA. Figure 3(c) 
indicates that no outflow emission is seen at positions more than ~ 11" from L1014-IRS in the 
30-m data and so no bright large scale outflow is present. In fact, there appears to be residual 
emission (cloud plus outflow) detected with the 30-m at the position of L1014-IRS some of which 
could be associated with the SMA outflow. This simple analysis suggests that any missing flux is 
limited in velocity and extent, and is less than a factor of 3 (Fig. 3(c)). 

The outflow properties are calculated in the standard manner (e.g., Cabrit & Bertout 1990). 
We assume a value for the excitation temperature of 20 K, but values in the range 10-50 K modify 
the outflow calculations only slightly, less than a factor of two. Two methods are typically used to 
calculate upper and lower limits to the outflow mass. For the lower limit the emission is assumed 
to be optically thin, the mass is calculated at each velocity channel and position where outflow 
emission is observed, and then summed. In the case of the SMA observations we are also missing 
some emission that has been filtered by the interferometer sampling, and so these values are strict 
lower limits. However, Fig. 3 suggests that the SMA observations do not suffer greatly from this 
problem. For the upper limits, a correction for line opacity is required, but we cannot estimate the 
opacity from our observations as we do not detect 13 CO J = 2^ 1, and the sensitivity of the line 
observations is insufficient to place any meaningful limits on the ratio of 13 CO to 12 CO emission. 
Thus, we derive only a lower limit to the outflow mass of Mfi ow ~ 1.4 x 10" 5 M , but note that, 
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if we correct for optical depth using typical values of 2-5 (Levreault 1988), then the upper limits 
to outflow mass and related properties will be correspondingly greater than our lower limits by 
similar factors. The mass estimate is consistent with the upper limit of 2 x 10" 3 M Q estimated by 
Crapsi et al. (2005) from their single-dish data. 

For properties relying on a knowledge of the outflow velocity (Momentum P, Energy E, 
Mechanical Luminosity L m , and Force F obs ), lower limits are determined by multiplying the mass 
in each velocity channel by the relative flow velocity (to the appropriate power) of that channel 
Vfiow =1 Vfiow-Visr I where Vfi ow is the observed flow velocity and V hr is the systemic velocity. Upper 
limits are found by assuming that the outflowing gas is moving at the maximum observed velocity 
V m , rather than the relative flow velocity of the channel, corrected for the outflow inclination. 
From modeling of the infrared scattered light nebula Huard et al. (2005) determine a semi-opening 
angle of 9 > 50°, and an inclination angle i > 60°. The correction factor 1 / cos(i-9) for inclination 
effects is then negligible (for i < 80°). The overlap between the blue and red lobes seen in Fig. 1(b) 
suggests that i is closer to 60° than > 75° (Cabrit & Bertout (1990)). Regardless of the true value 
for i, if 20° < i < 80° then the correction for inclination is negligible due to the large observed 
opening angle. Therefore, for the upper limit values we use the observed maximum flow velocity 
of 3.6 kms" 1 for V m , assuming no correction is needed for inclination. In this way we find P = 
2.4-5.0 x 10" 5 M Q kms" 1 ; E = 0.24-1.6 x 10" 4 M Q km 2 s" 2 ; L m = 0.34-2.8 x 10" 5 L Q ; and 
F obs = 3.4-7.1 x 10" 8 M Q kms" 1 yr" 1 . If correction factors for opacity of ~ 3 and missing flux of 
~ 3 (Fig. 3) are assumed, then the maximum values for P,E,L m and F obs will increase by about an 
order of magnitude. 

The outflow mass loss rate M out can be estimated directly from the mass and age tf [ow . Typi- 
cally tfi ow ~ td (the dynamical time) is assumed, which may not be reasonable (Parker et al. 1991); 
t d is a lower limit to the outflow age. Assuming tf tow ~ t d = 700 yr (assuming a lobe size of 540 
AU, see §4), and allowing as above for an increase in the mass by ~ 10 due to opacity and missing 
flux corrections, we find an upper limit of M out ~2x 10" 7 M Q yr" 1 . If instead tf tow ~ 10^ with no 
correction applied to the mass estimate we obtain M out ~2x 10" 9 M Q yr" 1 . 

4.2. Discussion 

The values calculated above show that the outflow is of low mass, and correspondingly weak 
in its momentum, energy, mechanical luminosity, and force, when compared to other outflows 
(Bontemps et al. 1996; Wu et al. 2004). In fact, L1014-IRS is the lowest luminosity source to 
date with a detectable bipolar molecular outflow, of 292 outflows with well determined bolometric 
luminosity L bo i (Wu et al. 2004). Our data suggests that L1014-IRS is not edge on, and so its L bo i 
estimate is reasonable. Even allowing for an upward revision of the outflow parameters due to 
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uncertainties in the CO opacity and missing flux, the outflow mass is at least an order of magnitude 
lower than any other outflow; the lowest previously reported has M ~ 10" 3 M Q (Bl-IRS; Wu et al. 
2004). The derived values of L bo[ , M, L m , and F obs are consistent with extrapolation of the relations 
among these quantities found for other outflows (e.g., Wu et al. 2004). L1014-IRS appears to be a 
weaker version of a typical outflow. The low values might suggest that the emission is in fact due 
to bound motions and not outflow emission, so using a size of 540 AU and velocity of 3.6 kms" 1 
we estimate an interior mass of ~4 M Q is required for bound motions, significantly larger than the 
core mass within the same radius (<0.5 M Q ; Huard et al. 2005). 

For comparison, the outflow from the low luminosity Class protostar IRAM 04191+1522 
(L ~ O.15L cf. O.O9L for L1014-IRS; Andre, Motte & Bacmann 1999) is large (> 14000 AU) 
with M > 0.03 M (Lee et al. 2002 from BIMA data) and F obs ~ 2 x 10 5 M Q kms" 1 yr" 1 (Andre 
et al. 1999). Unlike L1014-IRS, the CO outflow from IRAM 04191 is well detected with the IRAM 
30-m over a larger velocity range, up to 10 kms" 1 from the line center. The implied accretion 
rate, which follows directly from F obs , is over two orders of magnitude greater for IRAM 04191. 
In addition, IRAM 04191 shows clear evidence for infall in mm transitions of CS and is bright 
and extended in mm molecular lines such as N 2 H + (Belloche et al. 2002; Crapsi et al. 2005). 
So although L1014-IRS and IRAM 04191 have similar low luminosities and are embedded in 
dense cores, the properties of their outflows and cores show significant differences, which may 
signify different pathways to formation, and different end results. Its outflow properties suggest 
that L1014-IRS has a low accretion rate for an embedded source, and we can speculate that it is 
either (i) close to finishing its main accretion phase, or (ii) its accretion rate is intrinsically low. 
Either scenario implies that L1014-IRS will remain a very low mass object, and it may be the first 
example of an embedded proto-brown dwarf. Further information, such as a near-infrared spectral 
classification, is needed to address its evolutionary state and future. 

If L1014-IRS is indeed a proto-brown-dwarf, then our results suggest that BDs can form in a 
broadly similar manner to low mass stars (cf., the recent detection of outflowing gas from the young 
brown dwarf p Oph 102 by Whelan et al. (2005)), but that there might be significant differences in 
the details, by comparison with IRAM 04191. Alternatively the differences between L1014-IRS 
and IRAM 04191 might simply reflect differences in the mass and evolutionary state from one 
low-luminosity object to the next. 



5. Summary 

We have used the SMA to search for CO / = 2 — >• 1 outflow emission from L1014-IRS. The 
results are as follows: 
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1. A low velocity (< 4 kms ) compact bipolar CO outflow has been discovered centered on 
L1014-IRS, which directly associates L1014-IRS with the globule L1014 at ~200 pc and 
confirms its low luminosity (~0.09 L Q ). 

2. L1014-IRS is the lowest luminosity and perhaps the lowest mass source known to be driving 
a bipolar molecular outflow. 

3. The CO outflow is one of the smallest known in size (~500 AU), mass (< 10" 4 M ), and 
energetics (e.g., F obs < 10" 7 M Q kms" 1 yr" 1 ). These values are consistent with the trends 
observed between them for large samples. 
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Fig. 1. — Integrated intensity maps of CO 7 = 2^1 emission toward L1014, overlaid on a H 
band image (1.6 /im) from Huard et al. 2005. (a): Blue (red) contours represent blue (red) shifted 
emission between velocities (LSR) of 2.0 and 3.8 kms" 1 (4.9 and 7.0 kms -1 ). The images were 
made with natural weighting, and have been tapered with a 3" Gaussian. The contours are 2, 3, 4, 
... x the rms of 0.16 Jy/beam (0.25 Jy/beam) for the blue (red) shifted emission. The dotted lines 
are drawn from L1014-IRS to the peaks of the blue and red shifted emission that are spatially offset 
from L1014-IRS. The beam of 3"19 x 2."84, P.A. 24.°7 is shown at lower right, (b): As for (a), 
but with integration over 2.2 and 3.4 kms" 1 (blue) and 5.2 and 7.5 kms _1 (red). The images were 
made with Briggs' robust parameter of 0. The contours are 2, 3, 4, ... x the rms of 0.08 Jy/beam 
(0.12 Jy/beam) for the blue (red) shifted emission. The position and uncertainty of L1014-IRS is 
indicated with the white cross, and the beam of 1"17 x 0"97, RA. 88.° 1 is shown at lower left. 
The position-velocity cut shown in Fig. 2 which passes through L1014-IRS is indicated by the two 
short black lines either side of the outflow with RA. ~ 30°. 
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Fig. 2. — Position- Velocity (P-V) diagram for CO / = 2 — > 1 emission. The direction of the cut is 
indicated on Fig. 1(b), and offsets are relative to L1014-IRS. The cloud systemic velocity marked 
by the thick dashed line. The contours are 2, 3, 4, ... x the rms of 0.08 Jy/beam. 
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Fig. 3.— Comparison of the IRAM 30-m and SMA CO / = 2 — > 1 observations of L1014. Panel (a) 
shows the 30-m (11" beam) and SMA spectra (convolved with a 4" beam and scaled by the ratio of 
the 30-m and SMA beam areas) at the position of L1014-IRS. Panel (b) shows an enlarged version 
of a portion of panel (a) for clarity. Panel (c) shows again the SMA spectrum, and the difference 
("Residual") between the 30-m center spectrum and a mean spectrum made by combining spectra 
at 1 1" offset positions in a ring around the center. 



